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ABSTRACT 

The Nal Di line in the solar spectrum is sometimes attributed to the solar chromosphere. We 
study its formation in quiet-Sun network and internetwork. We first present high-resolution profile- 
resolved images taken in this line with the imaging spectrometer IBIS at the Dunn Solar Telescope 
and compare these to simultaneous chromospheric images taken in Call 8542 A and Ha. We then 
model Nal Di formation by performing 3D NLTE profile synthesis for a snapshot from a 3D radiation- 
magnetohydrodynamics simulation. We find that most Nal Di brightness is not chromospheric but 
samples the magnetic concentrations that make up the quiet-Sun network in the photosphere, well 
below the height where they merge into chromospheric canopies, with aureoles from 3D resonance 
scattering. The line core is sensitive to magneto- acoustic shocks in and near magnetic concentrations, 
where shocks occur deeper than elsewhere, and may provide evidence of heating sited deep within 
magnetic concentrations. 

Subject headings: Sun: photosphere - Sun: chromosphere - Sun: magnetic fields - Sun: faculae, plages 



1. INTRODUCTION 

The past decade saw a rapid advance of our un- 
derstanding of fine structure and dynamics in the 
lower atmosphere of the Sun by combination of high- 
resolution observations with increasingly realistic nu- 
merical simulations. Ex amples are the form ation 
of the solar granulati on dNordlund et al.l I2009D and 
rever s ed granulation flLeenaarts fc Wedemever-Bohml 
[20051 ICheung et al~ll2007h . and the mag netoconvective 
origin of intergranular bright points (e.g.. iShelvag et all 
|2004|) . Such studies focused mainly on photospheric 
structure, and assumed local thermodynamical equilib- 
rium (LTE) in both the numerical simulation and the 
subsequent spectral-line synthesis. 

The present frontier in this type of research is to 
study solar fine structure higher up in the solar atmo- 
sphere, including the chromosphere. There the mag- 
netic field fans out and becomes more homogeneous, but 
shocks cause fine-scale structure and dynamics; LTE is no 
longer valid. Observationally, high-resolution and rapid- 
cadence imaging is required in spectral lines that sample 
these higher layers, with sufficient spectral resolution to 
chart and follow their small-scale and fast profile varia- 
tions. At present only imaging Fabry-Perot spectrom- 
eters at ground-based telescopes with adaptive optics 
perform this task adequately, sampling the few chromo- 
spheric lines in the red part of the visible spectrum. 

In radiation-magnetohydrodynamics (MHD) simu- 
lations, such upward advance requires full three- 
dimensionality and tackling the challenging complexi- 
ties of non-LTE (NLTE) radiative transfer and non- 
equilibrium ionization/recombination balancing within 
the simulation. Post-processing in the form of subse- 
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quent synthesis of spectral diagnostics requires 3D NLTE 
radiative transfer, in some cases even with partial fre- 
quency distribution. In this paper we make such an up- 
ward move for the quiet-Sun atmospheric regime sampled 
by the the Nal Di resonance line at 5895.94 A. This line 
is often taken to be chromospheric, but we show here 
that it is much less so than Call 8542 A or Ha. 

Recently, observational comparisons using Fabry-Perot 
spectrometry were made between Ca II H & K and Ca II 
8542 A by iReardon et al.l (|2009f ) and between Call 
8542 A and Ha bv lCauzzi et al.l (|2009f ). iLeenaarts et all 
(|2009t ) compared Fabry-Perot observations in Call 
8542 A with numerical synthesis of this line from a 3D 
MHD simulation snapshot. The present paper is a com- 
panion to the latter study (hereinafter Paper I), using 
the same snapshot but turning to the Nal Di line. The 
major modeling advance in both papers is the incorpo- 
ration of 3D NLTE radiative transfer for the profile syn- 
thesis of the respective line. Complet e frequency redis- 
tribution is assumed and is in or der (jUitenbroekl 119891 : 
lUitenbroek fc Brulsl[T992T ). 

However, instantaneous LTE ionization equilibrium 
was assumed within the MHD simulation itself. Most 
importantly, it was assumed for hydrog en even though 
the studies bv iCarlsson fc Steinl (|2002h in ID and by 
ILeenaarts et al.l (|2007f ) in 2D have shown that this is 
not an accurate assumption. Compared to LTE, hydro- 
gen is underionized in hot shocks and overionized in cool 
post-shock conditions. Both are due to a large ioniza- 
tion timescale compared to the hydrodynamic evolution 
timescale. Improved treatment is highly desirable, but 
its implementation in 3D MHD simulation codes presents 
formidable challenges including severe processing-speed 
requirements. Paper I and the present paper therefore 
represent steps along the arduous obstacle course towards 
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realistic ab-initio simulation of the chromosphere, but 
not the final one. 

In solar spectrum atlases Nal Di is about five times 
weaker than Call 8542 A in equivalent width, but its 
line center reaches deeper, about 200 K lower in bright- 
ness temperature. Nal D x and Call 8542 A indeed dif- 
fer in various formation properties. Nal Di is a strongly 
scattering resonance line from an easily ionized minority 
stage; Call 8542 A is an excited line of the dominant cal- 
cium ionization stage and has larger source- function cou- 
pling to the local temperature. These differences make 
it worthwhile to employ and compare both lines. 

In the next section we present and discuss high- 
resolution spectrally-resolved imaging of a quiet-Sun area 
in Nal Di. We then use 3D NLTE synthesis of this line 
in the MHD simulation snapshot of Paper I to analyze 
its formation in quiet-Sun conditions (Section [3]). We 
discuss the results in Section [4] and conclude the paper 
in Section [5l 

2. OBSERVATIONS 
2.1. Data acquisition 

We use imaging spectroscopy data obtained with 
the Interferometric BIdimcnsional Spectrometer (IBIS) 
at the Dunn Solar Telescope (DST) of the U.S. Na- 
tional Solar Observatory/Sacramento Peak. IBIS con- 
sists of two Fabry-Perot interferometers mounted in 
series and delivers sequences of narrowband images 
with the spectral passband stepping quickly through 
multiple user-selectable sp ectral lines (jCavallinil 120061 ; 
iReardon fc Cav alliui 2008). In this particular sequence, 
obtained during 15:30-17:07 UT on April 18, 2008, IBIS 
sampled Nal Di with a passband of 23 mA FWHM, 
Call 8542 A with a passband of 43 mA FWHM, Ha 
with a passband of 22 mA FWHM, and Fe I 7090 A with 
a passband of 24 mA FWHM. Additional images were 
taken in the G band with an independent camera using 
an interference filter of 10 A F WHM. One of t he two 
adaptive-optics feeds at the DST (Rimmelc 2004]) served 
to stabilize the telescopic image and to correct wavefront 
deformations in real time; the G-band images were ad- 
ditionally improved through speckle reconstruction. The 
30 x 54 Mm 2 field of view covered a quiet solar area near 
the edge of an equatorial coronal hole at apparent disk 
center (heliocentric longitude 0.56° W, latitude 5.3° S) 
with a pixel scale of 0.083 arcscc/pixcl (close to critically 
sampling the 0.16 arcsec diffraction limit of the 76 cm 
telescope at Nal Di). As part of a service-mode observ- 
ing program with Hinodc a full-Stokes scan obtained with 
the Hinodc/SP spectropolarimeter in the Fel 6302.5 A 
line between 14:32-16:00 UT overlapped with a portion 
of the IBIS field of view. 

2.2. Results 

The IBIS images shown in Figures HHIl were taken be- 
tween 15:54:53-15:55:17 UT and were selected as the 
sharpest of the sequence. The narrowband spectral im- 
ages were additionally corrected by the application of de- 
stretching vectors derived from simultaneous broadband 
reference images. The dashed box in Figure [1] specifics 
the overlap with the Hinodc/SP Stokes scan shown in 
Figure [2j The solid box outlines the subfield selected for 



Figure [31 chosen because it covers a similar area as our 
simulation with a comparable amount of magnetic net- 
work, including the most prominent bright points mark- 
ing magnetic concentrations. However, only its right- 
hand half overlaps with the Hinode/SP scan. 

The four images in Figure [T] differ clearly in their por- 
trayal of the solar atmosphere, ranging from the pho- 
tospheric G-band image at left to the chromospheric Ha 
image at right. These two extremes display characteristic 
quiet-Sun scenes. The G-band image shows regular deep- 
photosphere granulation with relatively few bright points 
marking strong-field magnetic concentrations, such as 
the network patch within the solid box. The Ha image 
shows chromospheric fibrils (or alternatively "mottles", 
but we use the term "fibril" for any elongated Ha struc- 
ture), shorter and more irregularly shaped than would 
be the case in regions with larger magnetic activity. The 
fibrilar structure is more pronounced around the photo- 
spheric network. 

The two intermediate images differ significantly from 
both the outer ones and also between each other. The 
Call 8542 A image shows some of the same fibrils seen 
in Ha around denser network, but in the quiet inter- 
network it shows a dark background with small roundish 
bri ghtness patches tha t are mostly due to acoustic shocks 
(cf. lCauzzi elld1[2009l). T his regime was called "clapoti- 
sphere" bv lRuttenT i 19951 ). denoting sub-canopy domains 
with sufficiently weak field that acoustic and gravity 
waves can run up, become shocks, and interfere in in- 
tricate patterns. Similar p atterns occur in field-fr ee hy- 
drodynamical simulations ( Wedem cver et aLll2004f ). 

The Nal Di image (second panel) shows essentially no 
fibrilar structure, with only occasional glimpses of short 
and bright elongated structures. The brightest G-band 
bright points show up as co-spatial bright-point features, 
but with wide bright aureoles around them. The dark 
internetwork areas show vague pattern similarities with 
the Call 8542 A image. Thus, at first sight Nal D 1 
seems to sample the same internetwork waves, and pos- 
sibly clapotispheric shocks, with photospheric magnetic 
concentrations visible as bright points with aureoles. The 
sa me is shown in the similar but lower-resolution images 
of IA1 et all (120021 ). 

With the lack of fibrils and only a few strong bright 
points in this quiet-Sun region, the Nal Di image ap- 
pears rather featureless. This apparent blandness is a 
solar property, since the broadband images taken simul- 
taneously with the spectral images show excellent sharp- 
ness, especially near the adaptive-optics l ock point (lo- 
cated in the subfield shown in Figure [3]) . IRighini et alj 
(|2009f ) have shown that IBIS achieves near diffraction- 
limited image quality in its narrowband channel. Only 
small-scale features with very low contrast might drown 
in the photon noise set by the average line-core intensity 
of 2000 photons/pixel during the 50 ms exposures. The 
variation in the appearance of the bright aureoles among 
different bright points and at different wavelengths shows 
that they do not result from scattering in the instrument 
or the Earth's atmosphere. 

Figure [2] confirms the identification of Nal Di bright 
points as magnetic concentrations. This part of 
the Hinode/SP Stokes scan lasted from 15:51:23 until 
15:59:38 UT, so that it is not fully co-temporal with the 
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Fig. 1. — Images from the DST. Panels: (a) G band, (b) Nal Di, (c) Call 8542 A, (d) Ha. The latter three are constructed from the 
IBIS profile samplings by determining the intensity of the line profile minimum per pixel, in order to compensate for Dopplcrshifts from 
vertical bulk motions. The dashed box shows the overlap with the Hinode/SP magnetograms in Figure [2] The solid box is the subficld 
selected in Figure [3] 



Nal Di images. Nevertheless, there is remarkable cor- 
respondence between the brightest Nal T>i features and 
strong longitudinal field, much more than with the (gen- 
erally much weaker) transversal field. 

Figure[3]contains enlargements of the selected subfield. 
The magnetogram in panel (c) shows the signed longitu- 
dinal flux density but has only partial overlap. Panel 
(d) adds a filtergram display of the Nal Di intensity, 
i.e., with the IBIS passband at the nominal rest wave- 
length of the line, as compared to the image generated 
using local Dopplcr compensation in panel (c). The two 
panels have the same greyscale. The brightest features 
are brighter in panel (d) due to large central downflows 
which can be seen in panel (f) as three tiny very dark 
patches in the middle of magnetic concentrations. Three 
other magnetic concentrations display comparably nar- 
row updrafts. These motions Doppler-shift the line core 
out of the passband, causing a dditio nal brightening, as 
also pointed out by IA1 et all (|2002t ). Such mixing of 
Dopplershift and intensity modulation generally affects 
standard fixed-wavelength narrow-band filtergrams. We 
therefore plot Doppler-compensated line-center intensi- 
ties rather than rest-wavelength line-center intensities, 
so that brightness variations are more indicative of the 
source function than of the velocity field in the line form- 
ing region. The larger-scale patterning in panel (f) has 
no obvious correspondence to the other panels. 

The upshot of this section is that the NalDI core 
brightness seems to coincide more with the photospheric 
Stokes-U magnetogram signal than with fibrilar struc- 
ture seen in Call 8542 A and H a. This sug g estion 
agrees with the earlier conclusion of iCauzzi et al.l (|2000f ) 
from lower-resolution (spatial and spectral) data that the 
Nal D2 core intensity correlates well with the apparent 
magnetic flux density in network (above 50 Mx cm" 2 
in their Figure Al). Figure Q] indicates that Nal Di in 



the internetwork may show some clapotispheric intensity 
modulation caused by waves and shocks, mimicking those 
in Call 8542 A, but none of the chromospheric structure 
seen in Ha. 

3. SIMULATIONS 
3.1. MHD simulation snapshot 

We used a snapshot from a radiation-MHD simu- 
lation computed wit h the Oslo Stagger Code (OSC, 
IHansteen et al.ll2007D . The same snapshot was used in 
Paper I. OSC includes an PTE equation of state based on 
the Uppsala Opacity Package (|Gustafssonlll973[ ) and uses 
radiati ve transfer with multi-group opacity binning fol- 
lowing iNordlundl (|1982f ) and NLTE scattering following 
ISkartlienn 2000f) to compute radiative losses in the pho- 
tosphere and low chromosphere. Radiative losses in the 
middle and upper chromosphere are evaluated with an 
escape probability method including continua and lines 
of hydrogen plus lines of Call. It was calibrated with the 
ID RADYN code (e.g.. ICarlsson fc Steinl[l997l) . Higher 
up, OSC employs optically thin radiative cooling. Ther- 
mal conduction along magnetic field lines is tak e n into 
account. More detail is given in IHansteen et al.1 (|2007f ) 
and Martmez-Sykora et al. (2008, 2009). 

The simulation used here has a grid of 256 x 128 x 160, 
corresponding to a physical size of 16.6 x 8.3 x 15.5 Mm 3 . 
It extends from 1.5 Mm below the photosphere to 14 Mm 
above it, i.e., from the upper convection zone into the 
corona. The grid has 65 km horizontal resolution. The 
grid is non-uniform in the vertical direction; its resolution 
varies from 32 km in the convection zone to 440 km in 
the corona. The snapshot contains bipolar magnetic field 
with a mean strength of 150 gauss in the photosphere. 

3.2. Radiative transfer 

For the radiative transfer computation of the emergent 
Nal Di intensity the atmosphere was interpolated in the 



4 




If » H If 20 Z4U X Z4 
n[Hm] ElMftil n [Mm] 

Fig. 2. — Nal Di profile-minimum intensity compared with lon- 
gitudinal and transverse magnetograms for the overlap between 
the IBIS and Hinode/SP fields of view. Panel (a): the intensity 
of the Nal Di profile minimum per pixel, as in Figure [T] Panel 

(b) : apparent longitudinal flux density. Panel (c): apparent trans- 
verse flux density. The flux density greyscale ranges from -441 to 
+629 Mx cm -2 for panel (b) and to +338 Mx cm -2 for panel 

(c) . The average measured flux density over the field-of-view is 70 
Mx cm -2 . 

vertical direction onto a grid with 153 points covering a 
height range between —0.5 and 2 Mm relative to con- 
tinuum optical depth T5000 = 1, covering the formation 
range of all pertinent sodium transitions. The electron 
density was computed assuming LTE ionization for all 
relevant species, including hydrogen, consistent with the 
LTE equation-of-state used in OSC. T he angle quadra - 
ture was computed using the A4 set of lCarlsonl (|1963l ). 
employing 3 rays per octant. The gas motions in the 
simulation were taken into account without adding addi- 
tional microturbulence. 

The radiative t ransfer compu t ation was performed 
using MULTI3D (jBotnenl Il997t ILeenaarts fc Carlsson 
120091) a code based on the ID code MULTI of iCarlsson 
(|1986| ) and including the same physics, but with a 3D 
short-characteristics solver allowing evaluation of the 3D 
radiation field. MULTI3D has been MPI-parallcllized us- 
ing domain decomposition so that it runs efficiently on 
supercomputers. 

Since 3D radiative transfer remains very demanding, 
the actual Na I Grotrian diagram was simplifi ed to a min- 
imal m odel atom. Inspired by the results of iBruls et alj 
(fl99l . we constructed a model atom of only four levels, 



consisting of the ground state, the upper level of Na I Di , 
an extra level at 1 eV below the continuum, and the Nail 
ground state as continuum. We gave the extra level nor- 
mal collisional coupling to the two lower levels but enor- 
mous (10 5 in excess of normal values) collisional bound- 
free coupling to the continuum and large (Gaunt factor 
2.2) radiative coupling to the upper level of Nal Di. We 
added this e xtra level to imple ment the suction process 
identified by IBruls et al.l (|1992f ) which draws population 
from the ion reservoir into the NalD ground state. We 
adjusted its strong couplings such that it has similar suc- 
tion power as the complete full-atom upper term struc- 
ture, which we determined by comparing the resulting 
Nal Di level populations to those of the full 1 8-leve l 
model atom specified in Table 2 of IBruls et al.l (I1992D 
for the standard ID FALC atmosphere of iFontenla et al.l 
(|1993| ). The result is shown in Figure 0] The emer- 
gent profiles in the upper panel are nearly equal. The 
lower panel compares the NLTE population departure 
coefficients 6; = ?if ,LTE /np TE for the lower and upper 
Nal Di level between the two model atoms. They match 
very well, especially over the most important range z = 
0.4 — 0.6 Mm where the line core becomes optically thin. 
The steep increase of the b\ departure coefficient is due 
to the dominance of photoionization over collisional ion- 
ization, but sufficient suction is needed to match it with 
the small atom. Slightly larger differences occur deeper 
down, where both small-atom curves miss the slight over- 
population in the FALC temperature minimum due to 
lack of photon suction in weaker lines. However, the ra- 
tio 62/61 remains nearly equal between the two models so 
that the Nal Di source functions are nearly equal; only 
the optical-depth buildup below the formation height of 
the line core differs somewhat. This simple model atom is 
accurate enough for our purposes. We used it to compute 
emergent Nal Di profiles from the simulation snapshot, 
both treating each column as a plan-parallel atmosphere 
and in full 3D, assuming complete redistribution which 
is a g ood assumption for this line (|Uitenbroek fc Brulsl 
[19921) . 

3.3. Snapshot properties 

Figure [5] displays simulation results in the form of hor- 
izontal cuts through the snapshot cube at four different 
heights, specified at the top of each column. 

The first row shows how the magnetic field is organized 
in small-scale intcrgranular magnetic concentrations, es- 
pecially at lane vertices, that expand with height as in 
the classical fluted fiuxtube paradigm fe.g.. lSpruitJll977l ; 
Solanki 1987, 1993). At z = 0.0 Mm (upper-left panel) 
they do not appear as roundish fluxtubes but rather 
as extended features similar to the G-band "ribbons" 
and "flowers" in observations with comparable resolutio n 
(jBerger et al.l 120041: IRouppe van der Voort et al.l 120051) . 
The expansion with height is striking (note the scale re- 
duction to seven times lower amplitude at z = 0.9 Mm). 
The actual field is bipolar with a roughly even split be- 
tween positive and negative concentrations, as shown in 
the first panel (where the field is mostly near- vertical) . 

The temperature cut at z = 0.0 Mm in the second row 
shows granulation. At z = 0.3 Mm the brightness pat- 
tern has changed into reversed granulation; the strongest 
magnetic concentrations stand out by being hotter than 
their surroundings at equal geometrical height. Higher 



■5 





^ 1 
- 

b 


t 

♦ ^ * \ 


c 




ID 



19 

eJM*1 



3D 



15 



If 

1 1 MM I 



UP 



If 

iIMral 



» 



Fig. 3. — Details of the subficld marked by the solid box in Figure Q] Panels: (a) G band, (b) Call 8542 A, (c) longitudinal field, (d) 
Nal Di intensity at the nominal line-center wavelength, (e) Nal Di intensity of the profile minimum per pixel with the same greyscale as 
panel (d), (f) Dopplcrshift of the profile minimum per pixel, with the greyscale ranging from —2 km s — 1 (dark, downflow) to +2 km s — 1 
(bright, upflow). 

tion of magnetic concentrations. At 2 = 0.0 Mm these are 
superimposed on the higher-density intergranular lanes 
that appear as negative of the temperature pattern. At 
z = 0.3 Mm the tops of the largest granules are dens- 
est. The magnetic evacuation persists at all sampled 
heights, spreading with the field strength but maintain- 
ing co-spatiality with the photosphcric topology up to 
z = 0.6 Mm. In the z = 0.9 Mm panel the gas den- 
sity shows slender high-density features that likely mark 
shock interference. 

The electron density (fourth row) mimics the temper- 
ature closely at all heights because both the hydrogen 
ionization balance and the electron-donor (Si, Fe, Al, 
Mg, Ca, Na) ionization balances are computed in LTE. 
At z = 0.3 Mm most magnetic concentrations have low 
electron density due to low overall gas density, but the 
hottest show excess electron density from donor ioniza- 
tion. The N c /N-n ratio of electron to hydrogen density 
(not shown) reaches values as low as 10~ 6 in the coolest 
clouds in the higher layers, well below the 10 -4 rela- 
tive abundance of the electron donor elements, imply- 
ing that even these are largely neutral. The thin high- 
temperature shocks around the cool clouds have large 
electron densities. 

The vertical velocity pattern in the bottom row shows 
the familiar pattern of granular upflow and downflow in 
the photosphere. Its amplitude diminishes rapidly with 
height. Higher up, the most striking features are large 
downdrafts in most of the stronger magnetic concentra- 
tions, in particular the strongest one at (a;, y) = (8.2, 7.2). 
They survive up to z = 0.9 Mm. The upflow patterns re- 
semble the gas density patterns. 

3.4. Emergent Na I D± profiles 

Figure [5] compares the spatially- averaged Nal Di pro- 
file from the simulation with the spatially-averaged IBIS 
observation and with a solar-atlas profile. The agreement 
between the IBIS and atlas profiles is excellent. The sim- 
ulation deviates in producing a slightly deeper core and 
narrower inner wings, with more pronounced knees. This 
lack of core broadening is similar to, but smaller than for 
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Fig. 4. — Comparison of Nal Di line formation in the FALC 
atmosphere between the 18-level (black curves) and 4-level (grey 
dashed curves) model atoms. Upper panel: emergent line profiles 
at £t = 0.95. Lower panel: NLTE population departure coefficients 
for the lower (n = l) and upper (n = 2) level of the line. 



up, fairly large dark patches ("clouds") of very cool gas 
occur, with hot white borders indicative of shocks. They 
occur preferentially near magnetic polarity changes. 
The gas density (third row) shows the familiar evacua- 
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Fig. 5. — Properties of the simulation snapshot. All panels show xy-cuts through the simulation cube at the heights specified at the 
column tops, with 2 = Mm corresponding to spatially-averaged T5000 = 1- Each row displays the quantity specified with its units at the 
right-hand side, from top to bottom: magnetic field strength, gas temperature, mass density, electron density, vertical velocity (bright for 
upfiow). Each panel is greyscaled to maximum contrast, with the range specified along each panel top. The top-left panel displays the 
magnetic field strength multiplied with the sign of B z (upwards bright) to show the polarity distribution. 

the Call 8542 A line in Figure 1 of Paper I. The latter 
was attributed, at least partially, to a lack of spatial res- 
olution in the microturbulcncc-frcc simulation. Since we 
use the same snapshot, this deficiency also applies here. 
However, its effect is probably smaller because Nal Di is 
formed lower than Call 8542 A whereas the temperature 
and density in Figure [5] show an increase in fine-structure 
amplitude with height. 

Figure [7J displays properties of the Nal Di spectral 
profiles that emerge vertically from the top of the simu- 
lation volume, in the form of images comparable to the 
lower row in Figure [3J The wing image in the first panel 
serves as low-photosphere indicator, showing the mag- 
netic concentrations as bright points. The blue Nal Di 
wing was not included in our earlier survey of proxy- 
magnetometry diagnostics (jLeenaarts et al.ll2006[ ). but it 
is comparable to the other strong-line wings there. The 
magnetic concentrations gain contrast because of weaker 
damping wings, and, in the case of Nal Di also because 
of enhanced ionization from partial evacuation. 
The second panel of Figure [7J emulates a narrow-band 
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Fig. 6. — Comparison of spatially-averaged observed and simu- 
lated Nal Di line profiles. Grey curve: averaged disk-center inten- 
sity from the NSO Fourier Transf orm Spectrometer (FTS) atlas 
calibrated by Ncckcl & Labs (1984). Red curve and red plus sign: 
IBIS observation. Black curve: simulation. Both the FTS atlas 
and the simulation provide absolute intensities. The IBIS data 
were scaled to match the FTS profile in the line core. 



7 




& S D 13 

n[MmJ 

Fig. 7. — Nal Di synthesis with the 3D solver, except for panel 
(c). The intensity maps show the vertically emergent intensity from 
the simulation. From top to bottom: (a) intensity in the blue wing 
at AA = —0.27 A from nominal line center; (b) intensity at the 
nominal line center (rest wavelength); (c) intensity in the profile 
minimum per pixel, using the ID solver; (d) intensity in the profile 
minimum per pixel; (e) Dopplershift of the profile minimum per 
pixel (upflows bright); (f) r = l height for the profile minimum per 
pixel. The dotted line near y = 2.0 Mm specifies the cut location 
used in Figure [8] The three symbols in the bottom panel specify 
the locations for which Nal Di formation breakdown is given in 
Figure [9] 



filtcrgram with the passband at nominal line center. The 
strongest magnetic concentration appears very bright, 
mostly through downdraft Dopplershift as demonstrated 
in Figure [9] The rest of the scene vanishes by setting the 
grey scale to include this feature. 

Panels (c) and (d) compare ID plan-parallel pcr- 
column radiative transfer with full 3D radiative trans- 
fer, plotting the profile-minimum intensity per pixel in 
order to compensate for bulk-motion Dopplershifts. The 
greyscales are identical. The scenes are comparable but 
the magnetic concentrations reach larger contrast in the 
ID case. In 3D, radiation scatters out of these and illu- 
minates their surroundings, as shown below. 

Panel (e) maps the Dopplershift of the line profile min- 
imum per pixel. The greyish small-amplitude bright 
and dark patterning resembles the z = 0.9 Mm vertical- 
velocity sampling in the last panel of Figure[5] The larger 
magnetic-concentration downdrafts mimic those in the 
z = 0.6 Mm panel most closely. Narrow bright structures 
due to upward-propagating shocks appear elsewhere. 

Finely structured bright shocks and dark clouds ap- 
pear abundantly in panel (f), which displays the geo- 
metrical height where radiation at the wavelength of the 
profile minimum per pixel reaches optical depth unity. 
The range is very large, varying from z = 0.5 Mm in 
the cool clouds to z = 1.1 Mm in the shock fronts. 
This rugged height-of-formation pattern is effectively a 
map of the Nal/Nall ionization-recombination balance 
and is closely similar to the electron density pattern at 
z = 0.9 Mm in Figure El adding the slender high-density 
features in the panel above it. 

In the cool clouds sodium is nearly completely ionized. 
This counterintuitive result (one would expect sodium to 
be largely neutral at these temperatures and densities) is 
explained as follows: The electron density in the model 
atmosphere is computed in LTE and the cool clouds 
have thus a low electron density. The NLTE sodium 
ionization-recombination balance is set by the interplay 
of photoionization versus photorccombination requiring 
one electron. The photoionization rate per sodium atom 
is independent of temperature, as it is set by the radi- 
ation field. The recombination rate per sodium atom 
however depends on the electron density. In the low- 
temperature clouds the electron density is low, thus 
the recombination rate is low and as a consequence the 
ionization-recombination balance shifts towards Nail. 

The Nal Di r = 1 height map obtains its enormous 
variation and fine structure from the large variations in 
electron density affecting the sodium recombination. The 
Dopplershift map in panel (e) follows this r = 1 height 
variation; it samples the Dopplershift in magnetic con- 
centrations and cool clouds deeper in the atmosphere 
than in hot shocks.. 

Note that the r = 1 height variation would be even 
larger if we had not plotted r v = 1 for the Doppler- 
compensatcd profile minimum per pixel but t„ = 1 at 
the nominal line center as in a narrow-band filtcrgram. 
In that case Dopplershifts would cause even deeper line- 
center formation, as deep as z = 0.3 Mm in the down- 
drafts in the magnetic concentrations (Figure [9]). 

3.5. Comparison of 3D and ID modeling 

Figure [8] compares our 3D evaluation of Nal Di for- 
mation to the ID one by showing an xz-cut through the 
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Fig. 8. — Comparison between ID and 3D Nal Di line-center formation along the xz-cut near y = 2.0 Mm marked by the dotted line in 
Figure[7] Panels (c) to (h) have the t = height overplotted for the 3D (black) and the ID solver (grey). Panels: (a) emergent intensity for 
the 3D (black curve) and ID (grey curve) solver in radiation temperature units; (b) magnetic field strength; (c) contribution function for 
the 3D solver; (d) gas temperature; (e) ID radiation field; (f) 3D radiation field; (g) ID source function; (h) 3D source function; (i) ID 
source function on line-center T„ -scales per column; (j) 3D source function on line-center T„ -scales per column. 



simulation volume along the y = 2.0 Mm dotted line in 
Figure [7] where it passes through the second-brightest 
magnetic concentration in panels (b) and (c). The mag- 
netic concentration stands out in panel (b) of Figure [8j 
which also illustrates its upward field spreading, and in 
panel (d) as a relatively hot deep finger. The cut passes 
through a few more intergranular lanes with fairly strong 
field but no vertex concentration, and samples a high cool 
cloud at a; = 2.5 Mm and a low one at a; = 8.9 Mm. 

The first panel shows the vertically emergent inten- 
sity at the nominal line-center wavelength along the cut, 
given in equivalent brightness-temperature units to fa- 
cilitate comparison with panels (d)-(j). The 3D curve 
has an aureole around the magnetic concentration and 
is generally much brighter and smoother than the ID 
curve. Panel (c) displays the 3D vertical contribution 
function in color. The nominal line-center t„ = 1 heights 
are overplotted, also in panels (d)-(h) as reference; they 
do not differ much between 3D and ID. The line-center 
formation height varies from z 0.4 Mm in the magnetic 
concentration to z w 1.1 Mm in the hot rim of the deep 
cool cloud that appears as a filamentary feature in the 
z = 0.6 Mm temperature panel of Figure [5j The other 
upward spikes in panel (c) sample similar features, all 
with very wide contribution functions. The cool clouds 
cause r Va = 1 dips. 

The lower three rows of Figure [8] analyze Nal Di 



source function properties, at left for the ID case and 
at right for the 3D case. Panels (e) and (f) map the 
profile-averaged angle-averaged radiation field 

POO p27T P 7T 

J=— / I„(p(v-v )sm9 &9&<j)&v, (1) 

with I u the local specific intensity, tp(y — vq) the local 
Nal Di area-normalized extinction profile, v the fre- 
quency, vq the nominal line-center rest frequency, and 
9 and <fi specify the viewing angle in spherical coordi- 
nates. The other panels map the total source function 
at v = vq. It can be written, dropping frequency indices, 
as: 

S= (l-e)(l-r))J + e(l-r])B(T)+r){l-s)B* (2) 

where J represents the Nal Di photon reservoir, B{T) 
the thermal reservoir, and B* the reservoir for round- 
about paths (in particular recombination cascades) that 
may populate the upper level and hence produce line 
photons. B* is only a formal Planck function, not nec- 
essarily corresponding to the local temperature. The co- 
efficient e is the thermal line-photon destruction proba- 
bility by collisional deexcitation after photoexcitation in 
the Nal Di transition and through H photoionization 
as background continuum process. The profile weighting 
with the latter contribution makes the total source func- 
tion frequency-dependent, even though the line source 
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function is not for complete redistribution. The coeffi- 
cient 77 is the probability of line-photon conversion via 
other "interlocking" sodium transitions from the upper 
level (in particular photoionization) . 

Panels (e) and (f) are identical in deep layers where 
J « B. Higher up, the two panels appear identical also, 
in this case because their structure is dominated by the 
thermal broadening in the 99(1^ — 1^0) profile weighting. 
Towards the corona (note that the red color in panel 
(d) represents the cutoff at T = 10 4 K, the actual tem- 
peratures are higher) the line profile becomes so wide 
that J approaches the emergent-continuum value with 
^rari ~ 5100 K (see the 500 nm pa nel of Figure 36 o f 
IVernazza et~aT1 11981L or Figure 3 of iBruls et aT1ll992f ). 
In the ID case this part (orange) is patterned as ver- 
tical search-light beams above hot granules. They are 
washed out in the 3D case. Similar wash-out occurs in 
the blue-black regime with T ra d < 3600 K where Nal Di 
is formed; for example, the cooler areas above r = 1 in 
panel (d) become black in panel (e) but are filled with 
warmer radiation in panel (f). The scattering contribu- 
tion dominates the source function in panels (g) and (h) 
in the blue domain around 2 = 1.0 Mm, but higher up the 
e B(T) and rj B* contributions mix in, the first producing 
the thin bluish- white curve marking the temperature rise 
to the corona (where e is negligible), the second produc- 
ing the slanted greenish features below this curve. The 
e B(T) contribution is the same in ID and 3D because 
it is local, but there is marked difference between ID 
and 3D for the 77 B* contribution due to non-local ion- 
izing ultraviolet radiation. However, these contributions 
are too high up to affect the Nal Di brightness, as is 
demonstrated in panels (i) and (j) by plotting the source 
function per vertical column on the relevant part of the 
corresponding line-center t„ scale. The magnetic con- 
centration is the hottest feature crossing the t„ = 1 line; 

all other locations sample the blue J-dominated scat- 
tering regime, of which the columnar structure in ID 
explains the modulation of the corresponding intensity 
curve in panel (a). The 3D curve corresponds to the 
wash-out of this structure in panel (j). The 3D source 
function in panel (h) shows a slanted dynamic fibril orig- 
inating from the magnetic concentration, appearing as 
a substantial low-reaching green 77 B* contribution, but 
even this does not affect the last panel because t„ = 1 
deepens within the concentration. T hus, this 3D simu- 
lation confirms the ID conclusion of IBruls et all (|1992T ) 
that the source functions of the NalD lines are domi- 
nated by resonance scattering. 

3.6. Nal D\ formation in specific features 

The first quartet of Figure [9] is a graphical break- 
down of Nal Di formation for the brightest magnetic 
concentration marked with a plus sign in Figure [71 us- 
ing the format developed by Carlsson & Stein (1994, 
1997) for Call H2V grains. It rewrites the contribu- 
tion function to the vertically emergent intensity C7 = 
dI v /Az = exp(— t„) x„, with x„ the total extinction 
coefficient, into a product of the three terms XiV t Vj 
and t„ cxp(— r„) that are plotted as greyscale quanti- 
ties in the first three panels with Cj as their product 
in the fourth one. The three terms in this decomposi- 
tion demonstrate that a large contribution to intensity is 



made at those positions that have a large opacity at low 
optical depth, with a high source function around opti- 
cal depth unity. The line source function in the second 
panel is frequency-independent because of the assump- 
tion of complete redistribution. Strictly speaking one 
should plot the total source function here. However, to 
increase readability we avoid the cluttered modulation of 
the frequency-dependent total source function and plot 
the line source function instead. This is justified because 
the line source function becomes equal to the continuum 
source function in the deep-formed LTE wings, and dom- 
inates the total source function in the scattering line core. 

The emergent Nal Di profile in the last panel of the 
first quartet is very asymmetric. It shows a bright hump 
at nominal line center that causes the extraordinary 
brightening of this feature in the second panel of Fig- 
ure [71 The actual line core is redshifted over 7 km s . 
The formation mechanism of this profile is rather simi- 
lar to the formation of internetwork Call H2y grains. A 
large (10 km s _1 ) downdraft in high layers moves the core 
opacity redward. Its inward integration at the profile 
minimum reaches r = l at z = 0.5 Mm, where the source 
function is dominated by scattering and has uncoupled 
from the Planck function. The profile minimum is much 
darker than the emergent intensity at nominal line cen- 
ter. However, it still stands out in the Dopplcrshift- 
compensated panels of Figure [71 evidencing a relatively 
high source function notwithstanding the signature loss 
by scattering. Figure [Jj suggests that the same occurs in 
other magnetic concentrations. 

Several circumstances conspire to produce the bright 
hump in the profile at nominal line center. Because of 
the gradient in the downflow velocity there is a location 
with high Xvl T v at Az^ = and z = 0.3 Mm (upper- 
left panel). This location happens to lie close to optical 
depth unity at that frequency (lower-left panel). Most 
importantly, the source function is not yet completely 
set by scattering, and partially follows the increase of the 
Planck function above z = 0.1 Mm (dotted line, upper- 
right panel). This gives rise to the local maximum in 
the source function in the upper-right panel of the first 
quartet. All effects combined yield a high contribution 
function (lower-right panel) and thus a peak in the pro- 
file. There is no corresponding peak at the red side of 
the profile because Xvl T v (upper-left panel) is low where 
t v = 1 coincides with the location of the source function 
maximum because of the velocity gradient in the atmo- 
sphere. 

The second quartet of Figure [5] shows a similar break- 
down at the asterisk in the cool cloud in the bottom panel 
of Figure [7] In this case the profile asymmetry is minor, 
there are only a small downdraft below z — 0.3 Mm and 
a small updraft above this height. The Planck function 
(dotted) drops steeply to 2400 K at z = 1.0 Mm. The 
source function happens to follow it far out, but it ac- 
tually uncouples already at z = 0.2 Mm where scattering 
sets in. It levels out at a radiation temperature 3100 K 
near z = 0.6 Mm where t„ = 1. The resulting line core is 
very dark. 

The third quartet of Figure[9]shows another breakdown 
for the slender hot front marked with a rhombus in the 
bottom panel of Figure where Nal Di has maximum 
formation height with t„ = 1 at z = 1.1 Mm. The Planck 
function has an extended minimum for T 5000 K 
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Fig. 9. — Nal Di and Call 8542 A formation breakdown for the three locations marked in the bottom panel of Figure [7] The first three 
quartets are for Nal Di, the lower-right one for Call 8542 A. In each quartet, each greyscale image shows the quantity specified in its 
top- left corner as function of frequency from line center (in Dopplershift units) and simulation height z. Multiplication of the first three 
produces the intensity contribution function in the fourth panel of each quartet. A r„ = 1 curve (grey) and the vertical velocity (white) 
are overplotted in each panel, with a v z = line in the first panel for reference. The upper-right panel of each quartet also contains the 
Planck function (dotted) and the line source function (dashed), in equivalent temperature units specified along the top. The lower-right 
panel also contains the emergent intensity profile, as brightness temperature with the scale along the right-hand side. The first quartet is 
for the brightest magnetic concentration at the location marked with a plus in the bottom panel of Fig [7] The second quartet is for the 
cool cloud marked by the asterisk in the same figure. The third quartet is for the hot front marked by the rhombus. The final quartet is 
for the same location but shows the formation breakdown for Call 8542 A. 
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and then rises steeply with height, through 6000 K at 
z = l.l Mm, but this is not followed by the source func- 
tion which uncouples much deeper and has a radiation 
temperature of 3400 K at z = 1.1 Mm. Thus, although 
Nal Dx is formed very high at this location, it does not 
show the relatively high temperature there. Note the 
large difference in line-core formation height between the 
cool cloud of the second quartet and the hot front of third 
quartet. This is caused by the strong dependence of the 
Nal Di opacity on the electron density. The cool cloud 
has a low temperature, hence a low LTE electron density 
and thus a low opacity due to relative lack of photore- 
combination (see Sec l3.4[) . The hot front has a relatively 
high temperature throughout the whole column, so that 
the relatively high opacity results in a large formation 
height. 

For comparison, the fourth quartet shows the corre- 
sponding Call 8542 A formation breakdown at the loca- 
tion of the hot front, confirming that this line is formed 
higher and has better coupling to the temperature: this 
location is quite bright in Call 8542 A (Figures 3 and 4 
of Paper I). 

4. DISCUSSION 

4.1. LTE ionization versus non- equilibrium ionization 

The MHD simulation is quite realistic in the photo- 
sphere but arguably less so higher up. One obvious 
imperfection is the assumption of LTE Saha-Boltzmann 
partitioning for hydrogen and the electron-donor met- 
als. These are largely neutral in the simulation's cool 
clouds, but sodium gets largely ionized there in our sub- 
sequent NLTE statistical-equilibrium modeling for that 
element because its ionization is dominated by hot pho- 
tosphcric radiation. The same would be the case for 
the other electron don ors in such modeling. W orse, 
iCarlsson fc Steinl (|2002ft and lLeenaarts et all (|2007ft have 
shown that for hydrogen the assumption of statistical 
equilibrium fails also, requiring the evaluation of non- 
equilibrium ionization-recombination balancing in which 
the past history of a given gas parcel affects its state, in 
particular suffering slow recombination in the cool af- 
termath of a fast-ionizing hot shock passage. Similar 
slowness may affect other electron-donor elements, al- 
though hydrogen is likely the worst. Therefore, the low 
electron densities in the cool clouds in Figure [5j with 
their large effect on the Nal Di t = 1 height in the bot- 
tom panel of Figure [7J are likely artifacts. The actual 
electron density may remain as high as N /Nn « 10~ 2 
in shocked gas, causing up to 10 4 increase in H~ opac- 
ity and therefore increased radiative heating above the 
photosphere wherever the H bound-free continuum has 
J v > B v . So perhaps the cool clouds will not become 
so cool in the case of non-equilibrium ionization. On 
the other hand, slow balancing affects the energy bud- 
get intrinsically by not instantaneously converting ther- 
mal energy into ionization or releasing it in recombi- 
nation, resulting in larger tem perature contrast across 
shocks (ICarlsson fc Steinl I1992D and longer-lived cooler 
post-shock temperatures dLeenaarts et al.l I2007t) . It is 
hard to predict what cool clouds future simulations with 
improved ionization treatment will produce, but proba- 
bly these will have larger electron density than the cool 
clouds here, and possibly larger Nal Di opacity. 



4.2. Chromosphere or clapotisphere 

Since we use only a single snapshot of the simulation 
here, we cannot identify the nature of its fine structure 
in higher layers; that would require analysis of its tem- 
poral evoluti on. Nevertheless , comp arison with the 2D 
simulation of iLeenaarts et al 



tion of Martfncz-Sykora ct al 



and the 3D simula- 
suggests that many 
of the fine-scale features in Figures [SHU are due to shock 
dynamics, with slanted wave-guiding along field lines in 
the neighborhood of the magneti c concentrations pro- 
ducing short dynamic fibrils (see lHansteen et alj 120061 ; 
iDe Pontieu et all 120071: Paper I). The scene in the right- 
most column of Figure [5J consists of slender hot shock 
fronts and cool expansion clouds. There appear to be no 
long fibrils in the simulation that would be comparable 
to those typically observed near network in Ha and in 
Call 8542 A (Figure Q]). This is consistent with the ob- 
served Nal Di profile minimum intensity image shown 
in Fig. [T] which also shows almost no evidence of fibrils 
in this quiet region. 

The simulated Ca II 8542 A line-center image from the 
same snapshot in Figure 3 of Paper I does not contain 
such long fibrils either, and shows larger disparity with 
the observed line-center scene than occurs here between 
panel (c) of Figure [3] and panel (d) of Figurc[7j We there- 
fore conclude that the simulation does not contain the 
dense carpets of extended fibrils that constitute the mag- 
netically dominated chromosphere, but that it models 
what we call the clapotisphere quite well. The simulation 
is a fair description of what is observed in Nal Di but 
is deficient in Ha-type fibrils in Call 8542 A. Such defi- 
ciency doesn't really matter for Nal Di because its inten- 
sity is determined well below the clapotispheric regime 
anyhow: although its t„ = 1 height in the third quartet 
of Figure [9] is clapotispheric, the source function there is 
dominated by J photons originating much deeper. Only 
the core Dopplcrshift is set at clapotispheric heights, but 
only outside the transparent cool clouds. 

4.3. Magnetic concentrations in Nal Di and Call 
8542 A 

The magnetic concentrations appear bright in both the 
observations and the simulation, even after compensation 
for the Dopplershift imposed by the large downdrafts 
that tend to occur at their locations, also both in the 
observations and the simulation. In the simulation this 
brightness is due to higher-than-elsewherc temperatures 
in the upper photosphere (z = 0.3 Mm temperature panel 
of Figure O . 

Near and within magnetic concentrations, field-guided 
shocks tend to appear deeper than elsewhere, as 
wa s already sh o wn in the classical 2D simulation 
of iSteiner et alJ (|1998[ ) a nd confirmed l a ter by the 
dynamic fibril studies of lHansteen et al.1 (|2006f i and 
IDe Pontieu et all (|2007t) . The strong narrow downdrafts 
in most magnetic concentrations (bottom row of Fig- 
ure [5]) may be post-shock downdrafts. As diagnosed 
in the first quartet of Figure [H they cause magnetic- 
concentration brightening with a similar formation mech- 
anism as the asymmetric Ca II H & K violet- wing bright- 
ening in non-magnetic H2V and K2V grains. These 
sample the sub-canopy shocks in the internetwork at 
z 1 Mm. Nal Di is formed much lower and samples 
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only the deep-sited shocks in and near magnetic concen- 
trations. The formation of Call 8542 A is intermediate; 
the large brightness of the magnetic concentration in Fig- 
ure 4 of Paper I is similarly formed through higher-up 
downdraft. 

The breakdown for the magnetic concentration in Fig- 
ure [9] shows a raised temperature (Planck function in the 
second panel) over z = 0.3 — 0.8 Mm; it is the hottest lo- 
cation in the z~0.3 Mm temperature panel of Figure [5l 
The other magnetic concentrations have similarly raised 
temperatures at this height. The corresponding source 
function hump at z = 0.3 Mm, which produces the blue- 
wing peak at nominal line center, lies within the upper 
photosphere outside, in the reversed granulation regime 
(Figured]). 

The vertical velocity curve in the first quartet Figure 
(overplotted in all panels) displays a series of steps sug- 
gesting shock passages. If so, these are less well defined, 
with less obvio us response in temperatu re, than those in 
Figures 5-7 of iCarlsson fc Stem] (|1997l ) from their non- 
magnetic H2V grain simulation with the ID RADYN 
code. That code used an adaptive mesh and had much 
better vertical resolution than the fixed mesh of OSC. 
Also, in 3D shocks are likely not to run up straight 
even within magnetic concentrations. Analysis of the 
full time-dependent simulation will tell whether this ve- 
locity pattern indeed stems from shocks, and may estab- 
lish whether shock heating matters in setting this deep- 
seated temperature increase in magnetic concentrations, 
as suggested by the high correspondence between excess 
temperature and excess downdraft in magnetic concen- 
trations in the second column of Figure [5] 

4.4. Comparison between observation and simulation 

The observed and simulated scenes in panel (e) of Fig- 
ure [3] and panel (d) of Figure [7] show similar morphology. 
Within the simulation, Nal Di is primarily a magnetic- 
concentration mapper in its brightness distribution, sam- 
pling these in the upper photosphere and including wide 
scattered-light aureoles. Elsewhere, the line-center inten- 
sities are fully dominated by scattering, producing very 
dark source function sampling with much wash-out of 
thermal structure. Since these formation characteristics 
are primarily photospheric, where we believe our simu- 
lation to be realistic, we suggest that the same apply to 
the Sun. 

In the simulation only the core Dopplershifts, imposed 
at the last scattering, are likely to sample clapotispheric 
heights, and that only outside the cool clouds where 
sodium is fully ionized. The Dopplershift map in panel 
(f) of Figure[3]has similar slender magnetic-concentration 
downdrafts as in the companion panel (f) of Figure [7] 
but less other fine structure, in particular missing out on 
the considerable downflows in cool clouds and the sharp, 
slender upflow features from hot shocks that the simu- 
lation presents. These differences may be attributable 
to lack of resolution in the observations and/or overesti- 
mation of clapotispheric dynamics in the simulation. At 
this stage, we can't tell. 

5. CONCLUSIONS 

Based on the analysis of our simulation we conclude 
that the brightest quiet-Sun features observed in the core 
of Nal Di are also the deepest: magnetic concentrations 



sampled at z « 0.3 — 0.5 Mm. The lower height holds 
when one docs not compensate for Dopplershifts in the 
line-center sampling and a core shift causes additional 
brightening. This happens in particular when the mag- 
netic concentrations harbor higher-layer downdrafts pro- 
ducing bluc-wing peaks through the same line-formation 
mechanism that causes non-magnetic internetwork H2V 
and K2V grains in Ca II H & K. Such downdrafts occur in 
most magnetic concentrations in the simulation and seem 
to occur regularly also in the observations, although even 
our high-resolution IBIS imaging barely resolves them. 

The next brightest component in quiet-Sun Nal Di 
scenes are relatively bright aureoles surrounding the 
magnetic concentrations, due to resonance scattering of 
NalD radiation from the magnetic concentrations. This 
is very much a 3D phenomenon requiring 3D NLTE mod- 
eling, as demonstrated in Figure [8l 

Elsewhere, quiet-Sun Nal Di scenes are dark even 
though at certain locations the formation height may be 
located much higher up, around z = 1 Mm within hot 
clapotispheric shock fronts. Their high temperature has 
no response in the Nal Di source function due to the 
strong resonance scattering in this line, but it does for 
the Call 8542 A line which combines higher formation 
with better temperature sensitivity. 

The upshot is that in quiet-Sun areas the Nal Di 
core displays primarily an upper-photosphere rendering 
of magnetic concentrations on the solar surface, with 
roughly the same surface patterning as in unsigned mag- 
netograms or in proxy diagnostics such as G-band and 
Ha-wing bright points, but contaminated with bright- 
ening from downdrafts and 3D-scattered aureoles. This 
similarity explains the correlation between photospheric 
magnetic flux density and N alDi core brightness re- 
ported bv lCauzzi et al.l (|2000l ). and it makes Nal Di just 
another photospheric-field indicator in irradiance moni- 
toring. It is not a chromosphere diagnostic, not even a 
useful clapotisphcrc one. Only usage of information im- 
posed at the last scattering (Dopplershift, polarization) 
seems of interest, but that suffers from the large variation 
in the height of the r=l surface at line center. 

The magnetic concentrations in the simulation are hot- 
ter at the z = 0.3 — 0.5 Mm Nal Di line-center sampling 
height than the outside upper photosphere with reversed 
granulation. Further analysis of the full time-dependent 
simulation may clarify how this heating comes about; 
Nal Di may be useful to constrain the mechanism. 
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